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Abstract 

Nonlinear tunneling current through a quantum dot (an Anderson impurity 
system) subject to both constant and alternating electric fields is studied in 
the Kondo regime. A systematic diagram technique is developed for pertur- 
bation study of the current in physical systems out of equilibrium governed 
by time - dependent Hamiltonians of the Anderson and the Kondo models. 
The ensuing calculations prove to be too complicated for the Anderson model, 
and hence, a mapping on an effective Kondo problem is called for. This is 
achieved by constructing a time - dependent version of the Schrieffer - Wolff 
transformation. Perturbation expansion of the current is then carried out up 
to third order in the Kondo coupling J yielding a set of remarkably simple 
analytical expressions for the current. The zero - bias anomaly of the direct 
current (DC) differential conductance is shown to be suppressed by the alter- 
nating field while side peaks develop at finite source - drain voltage. Both the 
direct component and the first harmonics of the time - dependent response 
are equally enhanced due to the Kondo effect, while amplitudes of higher har- 
monics are shown to be relatively small. A "zero alternating bias anomaly" 
is found in the alternating current (AC) differential conductance, that is, it 
peaks around zero alternating bias. This peak is suppressed by the constant 



bias. No side peaks show up in the differential alternating - conductance but 
their counterpart is found in the derivative of the A C with respect to the di- 
rect bias. The results pertaining to nonlinear response are shown to be valid 
also below the Kondo temperature. 
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I. INTRODUCTION 



An attractive research direction in contemporary condensed matter physics seems to 
be the study of non-equilibrium many - body phenomena. A promising experimental and 
theoretical framework for investigating this topic is provided by the physics of quantum 
dot systems. The reason for that is clear, namely, quantum dots are fabricated and their 
properties can be elucidated by present day experimental techniques. Indeed, recent ex- 
periments on electron transport in quantum dots at low temperatures reveal signatures of 
coherent many - body physics such as the emergence of zero bias anomaly in current - volt- 
age characteristicsEHi which is due to the formation of a many body resonance. At the same 
time, the underlying theoretical models are of sufficient simplicity, so that one encounters a 
rare occasion where one has an experimentally accessible non-equilibrium quantum system 
that is also amenable to reliable and controllable theoretical approaches. 

So far, the main effort in the physics of quantum dots has been devoted to the study of 
different phenomena emerging from the presence of large Coulomb interaction in the dot. 
In this contexts, the most familiar and simplest topic is the Coulomb blockade. Its essence 
is simply encoded as a capacitance effect: every extra electron coming into the dot has to 
overcome a charging energy e 2 /C, where C is the capacitance of the dot. If the gate voltage 
is not tuned to supply this energy, the tunneling is (Coulomb) blocked. Yet, we believe that 
essentially genuine many - body aspects of resonant tunneling might be better revealed in 
the physics that goes substantially beyond the simplified Coulomb blockade picture. By this 
we mean the Kondo effect and other facets of strongly correlated electronic systems. They 
result from an intricate combination of electron - electron interaction and tunneling. 

Whereas the Kondo effect in bulk materials has been a thought - inspiring subject of 
research for more than three decades (for review, see Ref. its emergence in quantum 
dot physics appears to be relatively new. Yet, it proves to be equally thought - inspiring. 
In particular, it opens a road to explore the non- equilibrium Kondo physics. Its hallmark 
is the zero - bias anomaly, that is, an appearance of a large narrow peak in the differential 
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conductance around zero bias. Pertinent experiments have been carried out on crossed - 
wire tungsten junctionsi, quenched lithographic point contactsl'l, metal and metallic glass 
break junctionsi^ and, recently, quantum dotsEHl 

In parallel, progress has also been recorded in numerous theoretical works,0~@. Among 
other directions of research, interest is focused on time - dependent aspects of the Kondo 
physicf&i. Although an actual experimental research has not yet been carried out, 
measurement of alternating tunneling current in the appropriate range of frequencies proves 
to be feasibleHi,ii. Moreover, application of an external alternating electric field non- 
trivially affects the dc, an observable which can actually be measure We hope that 

the pertinent experiments can be carried out in the near future. 

It is useful to briefly mention the main features of the Kondo effect in bulk systems and 
in quantum dots. The Kondo effect was first revealed in the early 30 's as an enhancement in 
the resistivity of certain metals with decreasing temperature. Being a puzzle for thirty years 
it was then attributed by KondcS to interaction between conducting electrons and magnetic 
impurity atoms. Subsequent investigations proved that at low temperatures a hybrid state 
is formed consisting of the conducting electrons assembled around an impurity so that at 
zero temperature the magnetic moment of the impurity is completely screened. 

Within the realm of quantum dot physics, the impurity spin is represented by the spin 
of a single electron which is virtually locked in a deep level of the quantum dot. The 
formation of the hybrid state (either in bulk systems or in quantum dots) is accompanied 
by an appearance of a narrow peak in the interacting density of states of the impurity 
(or dot electron) close to the chemical potential of the band (lead) electrons. This peak 
is termed as an "Abrikosov - Suhl resonance" or "Kondo resonance" (see Ref. f|, pp.109, 
127 - 132, 210). For quantum dots in a static non-equilibrium situation, when two chemical 
potentials are present in the system (one for every lead), the pertinent peak in the interacting 
density of states splits, under certain conditions, into two peaks, one at each chemical 
potentiailllifflP (see Fig. |a). 

The main object of our study in the present work is the response of a quantum dot to the 



application of a strong time - dependent external electric field. For simplicity it is assumed 
that the field is monochromatic, whose frequency Q is in the range of tens of gigahertz. The 
relevant Physics implies an interesting extension of the underlying Kondo physics since, in a 
time - dependent field, every eigenstate of electrons in the leads is into a family of 

states whose (quasi) energies are separated by the photon energy hfl. Consequently, there 
is a family of Fermi seas, and, correspondingly, a family of chemical potentials in each lead. 
Thus, the Kondo peak in the density of states might split into a whole set of peaks (see Fig. 
|l]b). Numerical calculations of the indicate that this is indeed the case. We think that 
the occurrence of many separated Abrikosov-Suhl resonances can be tested experimentally. 
It will show up as side peaks in the (otherwise monotonic) I(V) curve displaying the zero 
bias anomaly of the direct current. Another interesting question is how the contributions 
of these peaks to the current, which are mutually coherent, interfere. This interference can 
show up only in the time - dependent currentS 

In this work, both direct and alternating tunneling currents are calculated within a 
perturbation scheme specially adapted for systems out of equilibrium. So far, calculations 
of the current through a Kondo system subject to a time - dependent bias were carried 
out using various assumptions and approximations. Here we carry out straightforward non- 
equilibrium perturbation expansion of the current in powers of the coupling constant between 
the quantum dot (or impurity atom) and the electrons in the conduction bands. Indeed, 
perturbation theory proved its usefulness in calculations related to the equilibrium Kondo 
model in bulk systems.il'11'0. Recently it has been used for quantum dots described by the 
Anderson model (in the Kondo regime) with a constant voltage biasEl. 

Let us then briefly list the main achievements of the present study. Our formulation 
starts from the time - dependent Anderson model which has already been shown to ade- 
quately describe the essential physics of a quantum dot (or a tunneling system based on 
an impurity atom) at voltages and frequencies less than the level spacing in the tunnel- 
ing region. We then develop the pertinent diagram perturbation technique which is able 
to treat non-equilibrium tunneling problems such as the time - dependent Anderson and 



Kondo models (for the former one it employs the slave - bosons method). A key point in 
the derivation is provided by combination of a specific approach suggested by ColemanEll 
to perform self consistent quantum and thermodynamical averaging in strongly correlated 
systems, with the Schwinger - Keldysh non-equilibrium Green functions formalism. As it 
turn out, calculations pertaining to the Anderson model appear to be rather cumbersome. 
At this stage we are therefore content by stating the rules for calculating diagrams and by 
pointing out a specific example where the Kondo behavior shows up. At the same time, 
we find it more practical to map the original time - dependent Anderson Hamiltonian on 
an effective time - dependent Kondo Hamiltonian. This is achieved by introducing a time 
dependent Schrieffer - Wolff transformation. Perturbation expansion of the current is then 
carried out within the Kondo model up to third order in the Kondo coupling J (sixth order 
in the tunneling coupling between conducting electrons and the dot). Remarkably simple 
analytical expressions are obtained for the whole spectrum of the tunneling current. It is 
shown that the zero - bias anomaly of the DC differential conductance is, in general, sup- 
pressed by an external alternating field, while side peaks develop at higher source drain 
bias. The nonlinear time - dependent current is found to be an interference sum of "direct - 
like" contributions, each one with an effective bias determined by the number of absorbed or 
emitted photons. The interference is shown to be rather destructive for all higher harmonics 
except the direct and the first ones. These two harmonics (and only them) are enhanced as 
a result of the Kondo effect while the other harmonics remain relatively small. In this sense 
the Kondo system behaves like a usual resistor (although the current is Kondo enhanced). 
Namely, direct and first harmonic input voltage lead to a direct and first harmonic current 
response. This result is demonstrated to be remarkably different from that for a non - in- 
teracting one - level system where all the harmonics emerge together. Moreover, it is shown 
to be valid also below the Kondo temperature. Finally, we also calculate the differential 
alternating conductance and find that like its direct analog it displays a zero alternating - 
bias anomaly. Inspecting the dependence of the AC on the direct bias reveals a non - trivial 
structure marked by side peaks. 
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The rest of the paper is organized as follows. In the next section, diagrammatic rules 
for the non-equilibrium perturbation expansion of the time - dependent Anderson model in 
the slave - boson representation are formulated. The diagrams are drawn up to sixth order 
in the tunneling coupling between the leads and the dot (Vkd) and some formal analytical 
expressions for the current are obtained. Further elaboration of the perturbation expansion 



in this model meets computational problems which today look severe. In section |I| a time - 
dependent version of the Schrieffer - Wolff transformation is developed in order to transform 
the time dependent Anderson model onto a Kondo - type model. In section |IV| perturbation 
expansion of the current (in the time - dependent Kondo model) is carried out in powers 
of the coupling strength (Jk'k)- An analytical expression for the current is obtained and 
the novel results arising from that expression are discussed. The paper is concluded with 
a summary which, in particular, includes some prospects for further research directions. 
In view of the inherent complexity of the pertinent formulation, we try to present it in a 
pedagogical style. 



II. PERTURBATION EXPANSION OF THE CURRENT IN THE 
TIME-DEPENDENT ANDERSON MODEL 

A. Bare Hamiltonian and parameter specification 

In this section we first introduce the time - dependent Anderson Hamiltonian and specify 
the range of parameters appropriate for the pertinent physical problem. Then we recall the 
slave - boson approach to the Anderson model and combine it with the Schwinger - Keldysh 
non-equilibrium Green function formalism in order to get an equation for the tunneling 
current. Then we develop a perturbation expansion of the current in powers of the tunneling 
strength (Vkd) up to sixth order in which the Kondo physics above the Kondo temperature 
(or in a strong enough external field) is unraveled. 

The time - dependent Anderson Hamiltonian takes the form, 
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H A = 



k£L,R;a 



(e k + A L(R) (t)) a\ a a Ka + ^ t d)(7 c\ a c d)(7 + \u + ( v kdol,<r c d,a 

k£L,R;a 



(1) 



Here a\ a {ak,a) creates (annihilates) an electron with momentum k and spin a in the left 
{k G L) or the right {k G R) lead, c d(T (c dj(T ) creates (annihilates) an electron with spin a in 
the dot, €k and e d ^ a are single - particle energies in the leads and the dot respectively, U is 
the Coulomb interaction energy in the dot and n a = c^Q )ff . The transfer matrix elements 
Vkd between the leads and the dot are assumed to be small compared with e dj(T and U . The 
external fields are included through potential shifts of the leads AL(R)(t), defined as, 



A L (R)(t) = <f>L(R) + W L (R) cos (ttt + a L(i?) ) . 



(2) 



The first term above describes a constant potential bias, while the second one is due to an 
alternating field, which, for simplicity, is assumed to be monochromatic. We note that the 
chemical potentials in the leads are shifted by the same amount A L ^(t) as the single-particle 
energies, hence the population of energy levels in the leads remains intact. 

As far as the value of the parameters is concerned, our attention here is focused on the 
Kondo regime which is determined by the conditions, 



td,a < 0, e di(T + U > 0, 



\Cd,<r\ , £d,a + U > T a , 



(3) 



where IV = 2tt EkeL, 



R 



Vkd fi(£d,a — £fc) are the widths of the energy levels in the dot. 
Furthermore, it is assumed that the external fields are not strong enough to draw the system 
out of this regime, so that: 



<f>L(R) < \£d,<r \ , ed,a + U and 



(4a) 



n,w L ,w R <\ + u. 



(4b) 



It should be stressed however that these conditions do not imply a linear response regime. 



The latter is defined by the conditions 



>L(R) 



, W L ( R ) < T while T < \e d ^\ , e d>(7 + U. 



At this point it is convenient to apply a canonical transformationoa on the Anderson 
model ([[]) whose purpose is to transfer the dependence on time into that part which contains 
a small parameter. The transformation is defined as follows: 



d 

H A = U' X H A U - ^- l —U (5) 
U(t) = exp {- J* fa [A L (h)N L + A R (h)N R ]} . 

where Nl(r) = J2a,keL(R) a k,a(t) a k,a(t)- The Hamiltonian Ha resulting from this transforma- 
tion reads, 

H A = E e kal,a a k,a + E e ^4, CT C d,a + ~U E n d,<r n d,cr> + 
keL,R;a & cr,a'^a 

+ E (Vkd(t)al a c d , a + h.c). (6) 

feeL,-R;cr 

where 

VfcdCt) = V kd exp 1 1 J^Jt x A L{R) ft) j . (7) 

Usually, in tunneling systems the barriers have low transparency. There- 
fore it is convenient to consider the tunneling part of the Hamiltonian, i.e. 
J2keL,R,a (ykd(t) a k,a c d,a + h.c.Y as a perturbationEll. It is well known, however, that if the 
interaction J2a,a'^a n d,and,a' is kept in the "unperturbed" Hamiltonian the Wick's theo- 
rem can not be applied. In order to circumvent this problem we assume that the energies 
associated with the direct and alternating voltages as well as the pertinent frequencies are 
smaller than the Coulomb interaction energy U in the dot. In other words, the assumption 
U — *■ oo should be an excellent approximation. We then apply the method of slave (auxil- 
iary) bosonsS^Bil using a certain version of it which is due to ColemanS Accordingly, the 
ordinary electron operators in the dot q i(J , c\ a , which transform a singly occupied state into 
an empty one and vise versa, are factored into a boson operator and a fermion operatorlll: 

Cd,a = tff a (8) 

<L = Ml 



The slave boson operator (b) creates (annihilates) an empty state, while the slave fermion 
operator f a (/£) annihilates (creates) a singly occupied state. In this representation the 
Hamiltonian (||)becomes, 

H S b= E e k a{ a a kja + J2^aPja+ E (V kd (t)a{^ f a + h.c.) . (9) 

keL,R;a & a,k£L,R 

The Coulomb interaction does not appear in the Hamiltonian any more. Indeed, when U 
is infinite it completely eliminates the possibility of double occupancy of the dot. This 
projection is accomplished by including the annihilation operator b in equations (§), which 
prevents creation of the doubly occupied state, and the constraint that the total number of 
slave bosons and slave fermions, 

Q A = h \ h+ Y J stu, (io) 

must be equal to unity. The operator Qa commutes with the Hamiltonian (Q) so that the 
"charge" Qa is not changed during the evolution of the system. Therefore the requirement 
Qa = 1 does not really constrain the dynamics of the system. Rather, it assures that the 
initial state does not contain a doubly occupied state as a component. In order to enforce 
the condition Qa = 1 Coleman0 introduced a Lagrange multiplier (chemical potential) —A 
and considered the Hamiltonian Hsb + ^Qa- The calculations are to be done with finite A 
and then, at the end, A —>■ oo. If the tunneling part of Hsb + ^Qa is chosen as perturbation 
the remaining "unperturbed" Hamiltonian is now quadratic in the creation and annihilation 
operators. The Wick's theorem can then be applied and diagrammatic expansion is feasible. 
We notice that the unperturbed part has become a simple time - independent free - particle 
Hamiltonian while the perturbation (tunneling part) contains both interaction and time - 
dependence. 



B. The tunneling current 



The tunneling current from the left (right) lead into the central region is defined as 
the product of electron charge (— e) and the rate of change in the number of electrons in 
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that lead. The latter is obtained by commutating the number - of - electrons (Heisenberg) 
operator Num with the Hamiltonian (|9|). This yields 

I L(R)( t ) = -4 dNL ^ {t) ) = -i E VUt)Tr FiA {p Fi 40)a{jtp(t)f (7 (t)} + c.c. 1 (11) 

\ ' a,k£L(R) 

where p f a(0) is the density matrix of the system at a certain fixed time which is taken here 
at t = 0. Generally, it should include all the changes that the system has undergone since the 
tunneling was switched onS We recall that the above equation is written in the Heisenberg 
representation. The angle brackets on the left hand side mean, of course, quantum average 
over the physical Hilbert space. At infinite U it is only a subspace F^ of the full Hilbert 
space for the slave - boson Hamiltonian (^) which is determined by the condition Qa = 1. It 
is explicitly manifested on the right - hand side of the above expression by the subscript F^ 
which restricts the trace to the physical subspace F^. In the following we calculate directly 
the current fllTD without prior elaboration on the dot Green functions as it was done in 
Refs. ll[ 13yi4yn| - [l8l^(^ Unlike the equation for the which expresses it in terms of the 
interacting density of states in the dot, an expression for the AC in terms of the dot Green 
functionsllil involves two Green functions (G < and G r ) and integration over real time. We 
find it easier to establish a perturbation expansion for the current itself (rather than for 
the Green function). To this end we rewrite equation (|IID in the interaction representation 
using the "grand - canonical" Hamiltonian Hsb + ^Qa- The result is, 

ll(R)(t) = -J E VUt)Tv Ff {p^(0)S(-oo, +oo)f [al >a (t + 0p(t)f a (t)S(+oo, -oo)]} + 

a,k£L{R) 

(12) 

where T is the usual time - ordering operator, S^+oo, — oo) = Te l $-<*> HTdt is the usual S 
- matrix and S(— oo, +oo) = S^+oo, — oo), while H T = J2keL,R-,cr (^4d(^) a i i(J c rf,o- + h- c ) is 
the tunneling part of the Hamiltonian Hsb + ^Qa which is chosen as a perturbation. The 
operators a\ a {t + 0), b\t) and f a (t) in the above equation now appear in the interaction 
representation (this is self evident as indicated by the presence of the S - matrix). We 
notice that it does not really matter whether one defines the interaction representation 



using the Hamiltonian Hsb + ^Q 'a or Hsb because Qa commutes with all parts of Hsb- The 
factor S(— oo, +00) in equation ([12]) prevents derivation of Wick's theorem and subsequent 
development of the Feynman diagrams technique. Indeed, the derivation of Wick's theorem 
is based on commutation of operators which produces pairings, namely "contractions" (see 
Ref. Commutation of two operators can not be worked out if one of them is subject 
to the time - ordering operator while the other one is not. The idea of a closed time 
path was introducedHSS in order to treat the troublesome factor S(— oo, +00) in non- 
equilibrium systems. The normal time branch is continued and turned back so that this 
operator becomes a factor within the scope of the time - ordering operator (T p ) on the 



closed time - path. Equation dl2|) , then, can be expressed as, 



iI( R )(t) = -j E v M (t)TT FiA {p FiA (o)fMA t ^( t -)f^) s p]} + c ^ ( 13 ) 



means 



Here S p = e l h HTdtp } where t p is the variable on the closed time - path and f 
integration over it. The subscript "— " on t signifies that the instant of time is considered on 
the normal (forward) time branch, while the subscript "+" is used for the backward oriented 
time branch (see Fig. |2|). The choice of the subscripts "+" and "— " for the time arguments 
of the operators a\ a (t+), ${t-) and f a (t-) is consistent in order to assure that their ordering 
by the operator T p is the same as in Eq. flllf) . The procedure of transformation from the 
Heisenberg representation (Eq. (|HD) to the closed time path representation (Eq. fll3|) ) is 
quite familiaiSOS We note, however, that it has usually been employed in the definition 
of non-equilibrium Green functions while here it is applied for the current operator. 

Yet, perturbation expansion of an expression for the current in the form (^) can not be 
carried out because of the constraint on traces in the subspace Ff 4 . In order to get rid of it 
we apply the Coleman method^ which is feasible here since the equation (|13|) appears as an 
ordinary statistical average of the operator O a = T p [a k a {t+)tf f a (t-)S p ) in the subspace 
F^. Let us choose the initial distribution of the system (before the tunneling was switched 
on) in the full Hilbert space to be an equilibrium one (grand - canonical ensemble with the 
"chemical potential" —A). Then the initial density matrix is, 
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p(-oo) = e-^ Ho+XQA ^/Z G (X), (14) 

where H = Hsb — is the unperturbed part of the Hamiltonian (|9|) and Zq(X) = 
Tr |e _/3 ^ 0+A< 3 A )| is the grand - canonical partition function. Now let us consider the expec- 
tation values of the operators OaQa and Qa in the full Hilbert space: 

(0 A Q A ) =Tt{ p (0)O a Q a }, 
(Qa)=Ti{ P (0)Qa}, 

where p(0) is the density matrix of the system in the full Hilbert space at the zeroth instant 
of time. When the tunneling is switched on the density matrix evolves in time, so that 
p(0) differs from p(— oo). However, due to the fact that the operator Qa commutes with 
the Hamiltonian it can still be factorized into separate blocks for each subspace Fq with 
different number Qa- Thus, 

oo 

(O a Qa) = E {p Fi (0)O A }Q A e-^Z QA /Z G (X), (15) 

Qa=0 

oo 

(Qa) = E QAe^ XQA Z QA /Z G (\), 



! A 



=0 



where Pfa(0) is the density matrix for the subspace Fq at the zeroth instant of time, while 
Zq a ee Tr F A |e _/3i/ °|. It is easy to see that in the limit A — ► oo the ratio of the two 
expressions written above becomes the expectation value of the operator Oa in the physical 
subspace Ff: 

(O a ) f a ee Tip a { Pf a(0)O a } = lim ^p^- (16) 

It is clear that the operator Oa has zero expectation value in the subspace F because there 
are neither slave fermions nor slave bosons there. Then Qa can be dropped out of the 
numerator in flTE|). This results in our final expression for the current, 



h(R)(t) = —r E v kd(t) hm ■ — -r + c.c. (17) 

In this equation the averages are taken in the full Hilbert space. It is amenable for standard 
field - theoretical calculation techniques. 
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C. Diagrammatic expansion of the current 



In the rest of this section we show how to work out a perturbation expansion for the 
current starting from equation (|17]) in powers of the tunneling part of the Hamiltonian In 
principle, a careful calculation of the density matrix p(0) is required for the non-equilibrium 
perturbation expansion. Otherwise, truncation of the full set of diagrams to a finite number 
might lead to completely wrong results. Indeed, truncation of a perturbation set implies 
that the effect of the perturbation is small, i.e. higher order terms of the perturbation series 
are negligible. It is not always ensured by the small value of the perturbation parameter. 
Consider, for example, an empty dot connected very weakly to a reservoir. The smaller 
the tunneling between them is, the longer it takes to fill the dot, but eventually the dot 
becomes full. It means that the occupation of the dot changes by a finite amount although 
the coupling is infinitesimally small. This "long - times" perturbation problem appears in 
the energy representation as an infrared divergenc e§. Fortunately, it has been explicitly 
shownll! that in the Kondo limit (determined by the conditions @ and (|4a|) ) the calculation 
of the Kondo - type contribution to the DC up to sixth order in V can be done using slave 
boson and slave fermion populations for a disconnected dot. In calculating the AC we also 
have to impose the condition (]4b"|) . The calculation of the denominator in equation (|17D is, 
then, very simple: 

(Q A ) =n b + 2n f A -=^ [l + 2e~^^] e~ p \ (18) 

where 

n b = 1/ (exp(/3A) - 1) and (19) 
rif = 1/ (exp[/5(e diCT - n + A)] + 1) , 

are Bose and Fermi functions for the slave particles. We assume here, for simplicity, that 
the number of spin degrees of freedom a is equal to two. 

In order to study the numerator of equation ( |TTD we start from its perturbation expansion 
in powers of (i.e. in powers of V^(t)). Then we transform every term of the expansion 
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from the closed time - path to the single time branch (see Ref. |57|j59| -|6lD . That implies 
expression of the Green functions as matrices in Keldysh space (see equations ( p0|) below). 
The next step is a rotation in this space (see Ref. |59| , |5"8D resulting in the so called@ "physical 



representation" for these matrices. Finally, time - translation invariance of the unperturbed 
Green functions allows us to apply Fourier transform and work in the energy representation. 

The diagrams of forth and sixth order in Vkd(t) are shown in Fig. |||. Recall that dis- 
connected diagrams need not be considered in the Schwinger - Keldysh formalism. Contri- 
bution of a disconnected diagram can be factorized into contributions of the connected and 
disconnected parts. Summation over all possible disconnected parts yields (T P S P ^ that is 
(S(— oo, — oo)) which is equal to unity. Notice, that diagram "h" is a crossed diagram. It 
is not included in the ever used non - crossing approximation. We now formulate the basic 
rules for drawing the diagrams and writing down the corresponding analytical expressions. 

1. In order to obtain a diagram of the m th order in V, draw a circle (it appears as a polygon 
in the figures) consisting of m alternating slave - boson and slave - fermion lines going 
in the same direction. The number m has to be even because every tunneling vertex 
contains only one slave - boson and one slave - fermion operator. The diagram must 
include only one such circle, otherwise it would produce contribution to the current of 
first or higher order in powers of exp(— f3X) (see Ref. fy]) that vanishes at A — > oo. We 
represent the diagrams by closed circles rather than by open lines as it is usually done 
since the external operators tft and f a in equation (|i7|) are taken at the same time t_. 
Connect the vertices by lead - electron lines. Remember that only the two types of 
vertices shown in Fig. ^ are allowed. 

2. Introduce the following 2x2 (Keldysh) matrices for every lead - electron, slave - 
fermion and slave - boson lines respectively: 




(20) 
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d{v) 



d a (u) ^ 



[d r (u) d%v) ) 



The indices "r" , "a" and "c" stand to denote retarded, advanced and correlation Green 
functions. They are defined explicitly as follows: 

1 



r(a) / 



e - e k + /i ± Z7 
^ = [1-2/(6*)] [^(e) - g a k M 
1 



(21) 



v - £d,a + /i - A ± Z7 

= [1 - 2 exp (-/3(6 djCT - /i)) exp(-/?A)] - £»] 

1 



V) 



v — A ± «7 
cf (z/) = [1 - 2exp(-/3A)] [<f(z/) - d°(i/)] 



where /(e^) = 1/ (exp[/3(efc — /i)] + 1) is the Fermi function for lead electrons, the fac- 
tors exp[— (3{e^cr — fj)] exp(— (3\) and exp(— (3\) are the limiting forms of the Fermi and 
the Bose functions ( fL9|) respectively at large A. We note that different representations 
can be used for the matrices (|20|) (see Ref. |58|-|61|). Here we employ the so called!! 
"physical representation" . 

3. Introduce the following tensors for every internal vertex (represented by a closed circle 
in Fig. 131). For the vertex drawn in Fig. |k: 



V kd J s (W ik) /n) e isa Wr]% 



(22) 



For the vertex drawn in Fig. ||b: 



v; d j q (w (k) /n) e-^«^. 



(23) 



where J s and J q are Bessel functions. The index "(fc)" reminds us that, despite the fact 
that W and a do not depend on k, they depend on the lead to which k belongs. The 
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factors J s ( q )(W(k)/Q) exp(±is(g)ct( fc )) originate from the time - dependence of V k d(f)- 
Indeed, substituting (|2|) into (|7|) we find 



V kd {t) = v kd e^^ t+iW ^ /nsin ( nt+a ^)' /3 ^ 

00 ■ r 1 

= v kd J2 J s (ww/n)en*w t+sUt+sa w- p w\ 



(24) 



The time - dependent factors appearing in this expression will surface when we impose 
energy conservation (see below), the phases (3{k) cancel, while the rest enters equations 
( p2"l) and (^). The physical meaning of s is the number of photons emitted when an 
electron goes from the dot to a lead, while q is the number of photons absorbed by an 
electron going from a lead to the dot. Both s and q can assume negative values. The 
tensor rj is given by the following expressions: 

72 



V2 



V 01 / 



r lf.i 



(25) 



V 10 / 



It appears as a result of the transformation leading from from time integration over 
past-ward going time branch to integration over normal time axis which requires a 
change of sign. The above form is consistent with the physical representation for the 
Keldysh matrices and tensors^. 

For the external vertex (represented by an open circle in Fig. |3|) containing the oper- 
ators a k a {t+), tf(t-) and / CT (i_) of equation (fl7|) write the factor J s (Wl/VL) exp(isaL) 
as for the internal vertices but do not insert the tensor 77. Instead, close around the 
product of matrices and tensors introduced above by the row 2~ 1 / 2 (1, 1) at the £ - 



matrix and the vectors 





and j. 




V 1 ) 




I 1 / 



at the matrices g and d respectively 



as we show in the example (equation (2S)) below. The external vertex differs from the 
internal ones because the time - variable for every external operator a\ a (t+), b\tJ) or 
f a (t-) is chosen on a certain (+ or — ) branch and these time branches are different 
for a\ a , and f a . In an internal vertex, the time - variables for all the operators are 
chosen on a single time branch following by summation over the two time branches. 
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4. Conserve spin in every vertex. As for energy, it should be conserved in every internal 
vertex taking into account emission (absorption) of s (q) photons and the different 
values assumed by the static potential energy <j)r k \ on the left and right lead. The 
energy in the external vertex is not conserved. The pertinent energy difference is 
equal to the frequency of the current. Therefore, introduce for this vertex the factor 
8n,(q-s+q 1 -s 1 +...), where n has the meaning of the total number of absorbed photons. 
Summation over n will be carried out in equation (|27|) below. 

5. Sum over energy, momentum and spin. Do not specify to which lead the momenta 
belong. It will be taken care of later on. 

6. Sum over numbers of photons q, s, qi, . . .. 

7. Multiply the result by the factor 

-J-(- ? r- 1 (i) 3m/2 (-l) P (2 7 r)- m/2 (26) 
HQa) 

where F is the number of closed electron - fermion loops. The first factor in the above 
expression comes from equation ( |T?D , the second one is implied by the m — 1 order of 
the expansion of the S - matrix. As for the next two factors, recall that an application 
of the Wick's theorem results in pairings of creation and annihilation operators which 
are, then, expressed through Green functions. The last factor simply results from the 
Fourier transform. 

8. At this stage take the limit A — > oo. 

9. Let us denote the expression obtained through the rules that have been listed by the 
symbol T^(n). Then the current is given by the following equation: 

lI(R)(t) = ~ E E {R* PM™)] «» (nnt - a w>(fcl) ,...) (27) 

71 n (k 1 ),(k 2 ),...=L,R 

+ Im [T A (n)] sin (nQt - a (fc)i(fcl)i ...) } , 
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where a 



(fc).(ki),- 



[s - q)a L + (si - ?i)a( fcl ) + (s 2 - g2)«(fc 2 ) + • • •, while £ (jfcl 



),(k 2 ),...=L,R 



means summation over the leads to which k\, k 2 , etc. belong. There is no summation 
over (k). The electron line coming out of the external vertex has momentum k be- 
longing only to the left (right) lead since we are calculating current through the left 
(right) barrier. The other momenta run through both leads. 



As an example we write below the expression which is obtained using the rules listed 
above for the diagram drawn in Fig. |3]e: 

1 1 



YaM = lim — — 

V ' A^oo K{Qa) 16 (27r) 3 



E 



6, 



n,(q+qi+q2—s—si—S2) 



(28) 



q,qi,q2,s,si,S2 



■ Jn 



Wl 



J, 



J, 



(kl) 
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n 



Wi 



(kl) 



J, 



Wi 



(*a) 



duii dui 2 dv\ dv 2 



J, 



Wi 



*2 



n 



/ „ Vkd Vk!d Vk 2 d 

ki£.L\k ,k 2 &L ,R\a ,a 2 

£ (i, i),r MvTMi 2 (^ - v - hk.) 

i,j,m,...,i 5 ,j 5 ,m 5 

^(u,! - ( gi - s 1 )Q) V rj 3 gi 3 j(uj 1 -v l -<j )L -{q + q 1 + q2 -s 1 - s 2 )Q) 



qi n)d m ^(u) v rj 2 



■ d m ^{v x + (q 2 - * a )n)i^*(u*)i$ 8 sg* (u* ~v x - {fc2) - q 2 Q) 

(-A (A 



\ 1 / 



V 1 / 



In the above equation the tensor product is to be expanded followed by integration over 
internal variables ui\, u 2 , V\ and v 2 is to be carried out. Using the Mathematica program we 
have carried it out for the above equation and for the analogous one given by the diagram 
drawn in Fig. |3|f. Whereas the calculation of the current could not be completed (see below) 
it was possible to inspect the emergence of a Kondo behavior through the appearance of 
a term J2 P /Jfl where f(p) = f(e p ) is the Fermi function for lead electrons of quantum 
number p. It is well known that this term is characteristic for the Kondo effect (above the 
Kondo temperature). Accordingly, the resulting expression for diagram ||f does not exhibit a 
Kondo - type behavior (i.e. InT - divergence in the linear response) while the one obtained 
for diagram and represented by equation (j28|) does. The equations in their final form are 
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very long and cumbersome and will not be shown here. Moreover, expressions appearing in 
the intermediate stages are unusually long and often overflowing the memory of a typical 
present - day workstation. We stress that this occurs only when computing the AC (calcu- 
lations appropriate for the DC are much simpler). These manipulations pertaining to AC 
appear to be especially cumbersome for the crossed diagram |3|h. We argue that, physically, 
the main source of complications results from the admixture of the Kondo resonance with 
usual resonant tunneling which is present in the Anderson model. Although the latter con- 
tribution is exponentially small within the relevant range of parameters, governed by terms 
like exp (ed, a /(kT)) and may frequently be neglected in the Kondo regime, its omission can 
not be assumed in this calculation scheme before arriving at the final stages. We are there- 
fore content by having introduced the systematic calculation scheme for the general time - 
dependent Anderson model, and exposing some diagrams that manifest a Kondo behavior. 
With rapid improvement of workstation capacities and analytic software programs the eval- 
uation of all diagrams of sixth order presented here can be completed according to the rules 
specified above. With this final note we now move on to complete the calculation of the 
current within an effective theory which still captures most of the pertinent subtle physics. 



III. TIME - DEPENDENT SCHRIEFFER - WOLFF TRANSFORMATION 

In this section we transform the time - dependent Anderson Hamiltonian to a time - 
dependent Kondo Hamiltonian. We use a procedure similar to the one proposed by Schrieffer 
and WolflS for the time - independent Anderson model but extend it to time - dependent 
problems. 

As we did in the previous section we start from the canonical transformation (|5|) of the 
Anderson Hamiltonian ([]]) whose purpose is to transfer the dependence on time into the 
perturbation part. It results in the Hamiltonian Ha (Eq. ([])). Then we further transform 
this Hamiltonian as, 

H' A = e s H A e- s -te s -^e- s , (29) 
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where the operator S is defined according to a specific prescription. In the spirit of Ref. 



44| we expand the right hand side of Eq. ( p9| ) in powers of S. We notice that S does not 
commute with ^ so this manipulation should be done with care. Then we require the 
operator 5* to satisfy the equation, 

(99 

i— + H T + [S,H o } = 0. (30) 

Here Ht = J2keL,R-,a (Vkd(t)al a c dt(T + h.c) is the tunneling part of the Hamiltonian Ha and 
Ho = J2keL,R;a e kal,a a k,a + E CT ^d,aC dta c d>a + \U E ff , ff '^ n d]a n di(T > is the rest of it. If S is 
proportional to V (which is indeed the case as we show at the end of this section) then 
solving Eq. (^) eliminates from the Hamiltonian H' A terms which are of first order in V, 
leaving only terms of higher orders which are naturally smaller. Upon collecting terms of 
zero and second order in V the Hamiltonian H' A can be written in the following form: 

H> A = H + l -[S,H T \ (31) 

Following Ref.0 we omit terms of 3d and higher orders in V. Our transformation differs from 
that of Ref.@ due to the presence of the term —ie s -j^.e~ s in Eq. ( [ZUD and, correspondingly, 
the term iff in Eq. (|3U|). 

Now we turn to the task of solving Eq. (|30"|), and look for a solution of the form, 

S = S x - S\. (32) 

The operator Si, then, should satisfy the equation 

+ H%*+[S 1 ,H o ] = 0, (33) 



■ ^>i , ,,,,„/ 



dt 

where H^ ut = J2k<=L,R;a Vkdityal a c d,a is the part of the Hamiltonian responsible for tunneling 
out of the dot. The operator S is anti - hermitian which assures the hermiticity of H' A . 
We are now looking for a solution of equation (|33| ) in the form, 

n d ,-<j , (1 - ni-o-) 



k£L,R;a 



£fc — ed,cr — U €k — e d ,cr 
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a i,*Cd,a- (34) 



The sum appearing in the square brackets is the inverse of the operator (k,a = £k — Cd,a — 
Uud-a- Schrieffer and Wolff found A k;a (t) = V kd , where V kd iSj of course, time - independent. 

Substituting expression (0) into equation fl33|) we find that A^ a (t) must satisfy the 
following equation: 



iA k (t)Ckl - A Ku (t) + V kd (t) = 0. 



(35) 



A proper solution of this equation is 



A k>a (t) = V kd £ J s (W ik) /a)e l ^ +s ^ t+lsa ^ 



(k,a 



(36) 



where V kd = V kd exp -i {w {k) /Q)j sina (fe ) 



+ Stt + Ck,a 

. Recall that the symbol (k) means "L" or "R" 



depending on whether k belongs to the left or to the right lead, while J s (-q) are Bessel's 
functions. The general solution of Eq. ([35|) contains also a term: oc e - *^^*, but the re- 
quirement that A kt(T (t) should be time - independent in the absence of external potentials 
enforces the prefactor to vanish. 

Employing now equations ([H]), (0), fl34|) and (|36|) we obtain the desired form of the 
Hamiltonian, in exact correspondence with Ref.0, 



H' A = H + ( w k'k{t) - l/2J k , k (t)n d - a ) 4 v a fejl 

k,k'eL,R;cr 

+ 1/2 J k'k(t)%r_ a a k ,aC dt<J C d - a 
k,k'GL,R;a 

~ 12 (W kk (t) - 1/2 J kk (t)n d ,-a)Cd ja c d ,a 

k£L,R;a 

_1 / 4 12 [ J k'k(t)al_ a al )a c dtCr c d _ a + h.c. 



(37) 



where 



Jk>k(t) = VjfdVfc d exp i {4>{ k ') - t 



(38) 



n 
i 



+ 



i(s' — s)Qt + i(s'a( k i) — sa^ k )) 
1 



+ (p(k') + s'fi - efc + <P(k) + sVt- e d!a 

1 1 



efc' + 0(fe') + s'fi - e di<J - £7 e fc + 0( fc ) + sfi - e d)(7 - C7 y 
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s',s=— 00 



1 1 

+ 



The form of this Hamiltonian is the same as that of Ref.@ but the matrix elements Jy k {t) 
and Wk'kit) are evidently distinct. 

In the Kondo regime, the important subspace F^ of the full Hilbert space is that one 
for which the dot is occupied by one electron. The last term in the Hamiltonian (j37D is 
not relevant for this subspace. The third term in equation fl37|) can be absorbed into the 
definitions of e d>a and U. Besides, we have n rfcr ?v_ CT = and J2a n d,a = 1 in the subspace 
Ff, so that the Coulomb interaction term (present in Ho) vanishes and the one - particle 
energy term for the dot becomes a c - number. The remaining terms represent the Kondo 
(also called "s - d" ) Hamiltonian plus a potential scattering term, 

H K = e k a k,a°-k,a + ^ ( W k'k(t) - l/2J k i k {t)n d ,- a ) 4', CT a ^ + 

k£L,R\<j k,k'eL,R;a 

+ 1/2 Jk'k{t)a\, _ a a k ^c ] d a c d - a . (39) 

k,k'£L,R;a 

Two comments are in order here: (i) Any procedure toward calculation of physical quantities 
should take into account the fact that, out of the full Hilbert space, the system is projected 
onto a subspace F^ in which the dot is occupied by one (and only one) electron, (ii) At 
this stage one might be tempted to express the electron creation - annihilation operators 
in the dot through spin operators, thus arriving at the familiar formBiS of the Kondo 
Hamiltonian. But then one would realize that the spin operators do not obey the usual 
commutation rules. In order to overcome this obstacle, fictitious (auxiliary) fermions might 
be introduced^. But this leads one back to equation (|39|). In other words, auxiliary fermions 
which are sometimes regarded as artificial particles introduced to represent spins are real 
electrons in the dot (impurity atom) subject to the constraint specified in (i). 

Calculation of the tunneling current starting from the Kondo Hamiltonian fl3TS| ) is possible 
for arbitrary field strengths and frequency provided the inequalities (§) are satisfied. Yet, 
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inspecting a typical experimental setuplira one may consider somewhat weaker external fields 
and lower frequencies, so that 



(f) L(R) , Q, W L , W R < \e di<7 \ , e d , a + U. 



(40) 



Expressions ( P5| ) for Jk'kif) and Wy k (t) then significantly simplify. Indeed, at small W/VL 
the Bessel function J S (W/Q) rapidly decreases with increasing s. At large W/fl it decays 
strongly once s exceeds W/VL. Therefore we can restrict sQ to be less than or of the order 
of max(f2, W), that is, sQ <C \e dta \ , e dia + U. In the formation of the Kondo resonance the 
most important states are those with energies je^l <C \^d,a\ > e d,a + U (see Ref. 44). Then 
\Ck,a\ ~ \^d,a\, if nd -a = 0, or \( k}(7 \ « e d;(T + U, if n d _ a = 1 (see definition after Eq. ||, and 
recall that n,i_ ff is discrete). Therefore the conditions (|40D assure that 



(41) 



(recall that refers to 0^ or <p R ). Therefore we can neglect the term + sVL in the de- 



nominator of the expression ([36|). Its right hand side is greatly simplified, and now becomes, 



A Ka {t) « V kd (t), 



(42) 



where Vk d {t) is defined by Eq. (0). Expressions for the matrix elements Jyk and Wykif) 
simplify as well. First, they become time - independent if k and k' belong to the same lead. 
Moreover, they do not depend on potential shifts of each lead separately but only on their 
difference, A^ — A^. For a monochromatic potential difference between the leads we then 
define 

A LR = A L -A R = (f) dc + W cos (Vlt + a) . (43) 

Then the matrix elements Jk>k{t) and W k >k{t) can be expressed in quite a simple form: 

i ft 



Jk'k{t) 



w k , k {t) 



J fc / fc exp ± \_Ji x A LR {t x ) , iik' G L, keR 



Ji 



k'k 



W k 'k 



if k',k G L or k',k G R 

, if k' G L, k e R 
if k',k G L or k',k G R 



(44) 
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where 



7 ~ V V* f 1 + 1 1 1 

Jk'k = Vk'dV kd | h 



W k , k = \v k>d V; d ( + ) . (45) 

^ \£k — Cd,a e k' — £d,a ) 

are time - independent. The matrix elements J k 'k(t) and W k i k (t) for k' 6 R, k E L satisfy 
the hermiticity relations J kk >(t) = J k i k {t) and W kk '(t) = W k , k (t). We note that at small 
external fields ( ^0[) the transformation (determined by Eqs. fl34|) and (]42])) and the form of 



the matrix elements J k ' k {t) and W k / k (t) are very similar to those of Ref.0 although they are 
still different due to the time - dependence of V kd (t). We further note that the Hamiltonian 
(|39|) with the coupling constants ( [44]) can alternatively be obtained by application of the 
canonical transformation (|5|) to the usual Kondo Hamiltonian with time - dependence added 
only to the leads (while the coupling of conduction electrons to the impurity (dot) remains 
time - independent). This is true, however, only at small external fields ([40D . At stronger 
fields equations (p8|) must be used for computing the matrix elements J k ' k (t) and W k i k (t). 

In order to have a more compact form of the Hamiltonian we hereafter imply the limit of 
infinite U. It is obvious from equations ([44]) and (^) that this choice does not lead to any 
qualitative changes. Indeed, it just slightly affects the values of J k i k . The main advantage 
of this choice is that at U = oo the equality, W k > k = 1/2 J k i k holds, which eliminates one 
more parameter. Furthermore, in the subspace the following identity holds, namely, 
1 — n d _ a = n d ^ a . Then H K can be expressed as 

Hk = e kak,a a k,a + (46) 

k£L,R;a 

+ V 2 J2 J k'k(t) [al,-a a k,aC d ^c d ^ a + al a a k)a c da c d ^ . 

k,k'£L,R;a 

There is only one coupling constant J k ' k (t) in this expression which is equal for both coupling 
terms. Equation (|4*6| ) constitutes our final form of the time - dependent Hamiltonian which 
we use in the next section to obtain the tunneling current. 

In concluding this section we would like to point out that the main idea of the Schrieffer - 
Wolff transformation is based on a projection of the system out of the full Hilbert space onto 
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the subspace F± for which the dot is occupied by one (and only one) electron. Unoccupied 
and doubly - occupied subspaces are forbidden. At first glance it looks as a small reduction 
of space dimension. Yet, as we show in the next section it greatly simplifies the calculation of 
the tunneling current for time - dependent problems. The main physical reason seems to be 
the fact that by fixing the number of electrons in the dot one separates the Kondo resonance 
from the usual resonant tunneling. The latter is exponentially small in the Kondo regime but 
it is formally present in the Anderson model. Admixture of two different physical processes 
within the same calculation scheme seems to be the main source of complication. Another 
advantage of the Kondo Hamiltonian (^) (or (f46|) ) is readily seen from the expressions 
(j39p, (|3|), ([44]) and (|4^). Namely, the number of independent parameters is significantly 
reduced. Indeed, instead of six parameters controlling the external fields in the Anderson 
model (<f>L(R), Wl(R), and the phase difference — we are left with only three: <p dc , 
W and Q. Instead of three independent internal parameters Vkd, £d,a and U there are only 
two important combinations Jyk and Wyk- Moreover, it becomes obvious that the limit of 
infinite U does not imply any qualitative changes in the results. It allows us to get rid of 
one more parameter. 



IV. PERTURBATION EXPANSION OF THE CURRENT IN THE KONDO 

MODEL 

A. Expression for the tunneling current 

In this section we define the tunneling current using the time - dependent Kondo Hamil- 
tonian fl46|) and then develop a non-equilibrium technique to expand it in powers of the 
coupling strength Jyk- Although the details of calculation are substantially distinct from 
those used in section ||, the basic algorithm is quite similar. Of course, we do not need to 
introduce slave particles here, because the unperturbed part is bilinear in creation - annihi- 
lation operators and there is no obstacle in carrying out perturbation expansion in powers of 
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the interaction which contains a small parameter J^k- The starting point is, again, defining 



the current using commutation of the Hamiltonian (|4q) with the number - of - particles 
operator. It yields, in the Heisenberg representation, 

TP IP 

I(t) = -j ({Hk,N l }) = - ([H K ,N R ]) = (47) 



% H Im{j^ fe (t)Tr F x [ Pf k (0)4, i _ CT (t)a jt)CT (t)4 i(J (t)c d _ <T (t) 



h 



+ 



+ Jjfc/ fc (t)Tr F K [p F x (0)4, CT (t)a fcjCT (t)4 i(7 (t)Q i(T (t)]} , 



where p f k(0) is the density matrix of the system at the zeroth instant of time (compare 
with Eq. fllTD ). As usual, the average is taken over the physical subspace of the full 
Hilbert space for the Kondo model, and the subscript on the right - hand side implies 
that operators and traces are performed within this subspace. Unlike the physical subspace 
for the Anderson model (section [TT]) which we referred to as Ff, the subspace Ff" is defined 
by the condition that Qk = 1, where 

Qi^E c W- (48) 

Definition of the physical subspace for the Anderson model, i.e. the condition Qa = 1 (see 
Eq. (|IUD) fixed the total number of slave particles to be equal to unity, thus, allowing both 
unoccupied and single - occupied states of the dot. The condition Qk = 1, on the other 
hand, enforces single occupation of the dot. 

In order to get rid of the constraint to the subspace F± we adapt the method proposed by 
ColemanS for the analogous problem in the Anderson model, to be used also in the Kondo 
model. As in section |TJ we combine it with the Schwinger - Keldysh non - equilibrium 
Green function technique. First, we introduce a grand - canonical Hamiltonian Hk + \Qk- 
The limit of infinite A is to be taken at the end of the calculation. Then we go to the 
interaction representation considering the exchange interaction as a perturbation and the 
rest of the Hamiltonian, i.e. Y,keL r-<t e k% a a k,a + ^Qk as the unperturbed part. We notice 
that Qk can be freely added to or subtracted from the Hamiltonian in the definition of the 
interaction representation since it commutes with all parts of Hk, although it is important in 
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the statistical average as we proceed to show below. Within the interaction representation, 
equation (pE?D for the current now reads, 



->t+u n kleL k( - R . . 

+ J k 'k(t)Tr F K [p F A'(0)f p (4, (T (t + )afc i(T (t / f )4,a( t -)c ( i,<T(t-)5' P )]} , 



where S p and T p are, respectively, the S - matrix and the time - ordering operator on the 
closed time - path. In this equation and hereafter, the operators a\, a (t), cifc j0 .(t), c\ a {t) and 
Cd,a{t) are defined in the interaction representation. The choice of points t + , t_, t' + and 
t'_ on the closed time path (see Fig. |^) assures proper ordering of these operators by the 
operator T v . As in section [il] the initial distribution of the system in the full Hilbert space 
(before tunneling was switched on) is chosen to be an equilibrium one. It corresponds to a 
grand - canonical ensemble with the "chemical potential" —A. Instead of the operator Oa 
we consider here the operator Ok which is defined by the equation, 



+ 



i v<7 (t+H At'+HA^) ^-) + 4v(*+)^4)4,,(OM*-)) s P ] . m 



K = T n 



Considering expectation values of the operators OkQk and Qk within the full Hilbert space 



and repeating the steps leading from Eq. (|T3|) to Eq. ([17]) (see equations (|T5|) ~~ (0) and 
explanations therein) we obtain the following expression for the current: 



lit) = — — Im lim lim { 

4' , -a (t + ) a k ,a (f+ ) cL (t'_ ) C di - a (t_ ) S p 



(50) 



Jk'kit) 



T 



+ Jk'kv) 



T 



(Qk) 

a l>A t AakA t '+) c dA t 'A c dA t ASp 



+ 



(Qk) 



Since the averages in this equation are taken in the full Hilbert space it is amenable for 
perturbation expansion. 

As in section [Tl] and Ref. |TB, we assume that calculation of the tunneling current to 
lowest order in the tunneling strength which encodes the Kondo effect (i.e. J|, fc ) can be 
done using unperturbed populations of the energy levels. This assumption seems more 
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natural here because the number of fermions in the dot is completely fixed by the constraint 
to the subspace . It is not the case in the Anderson model where slave fermions can 
convert into slave bosons and vise versa within the same subspace F^. Moreover, two levels 
with different spins have equal occupation numbers after performing an ensemble (thermal) 
averaging. Therefore, in the Kondo model this assumption concerns only occupation of 
different energy levels in the reservoirs. The denominator in equation ( |50|) can therefore be 
easily calculated, that is, 

(Qk) = 2n d ^ 2e^\ (51) 

where n<j = l/[exp(/3A) + 1] is the Fermi function for the dot electrons in the grand - 
canonical ensemble. In fact, as far as the dot electrons are considered in the grand - canonical 
ensemble, they are not real particles any more. For example, one might notice that in the 
physical limit A — > oo their Fermi function tends to zero. They might be better called 
"auxiliary fermions" or "slave fermions". To avoid confusion with slave fermions of section 
il] we prefer to refer to them hereafter as "dot fermions" . 



B. Diagrammatic expansion of the current 

In the following we formulate a diagrammatic technique to expand the numerator of Eq. 
(|50|) and obtain the current. We skip the detailed derivation and present only its main steps 
and then present an explicit formulation of the pertinent diagrammatic rules. The strategy 



is to start from perturbation expansion of Eq. floTf ) on the closed time - path. Next, trans- 
formation of the resulting expression is performed, first to a single - time branch and, then, 
to the physical representations of the non-equilibrium Green functions. Time - translation 
invariance of the unperturbed Green functions allows us to use the Fourier transform (for 
every Green function separately) and write integrals in the energy representation. The rest 
of this procedure will be explained along the course of formulation below. 

Equation ( |5T| ) expresses the current in terms of a pair of two - particle Green functions. 
Therefore, in this section we prefer to draw the diagrams in the standard form excepted for 
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Green functions in textbooks (see Fig. |6|). Alternatively, we could connect all the external 
lines in one "external vertex" as we did in section m thus, obtaining diagrams of a circular 
shape. The difference is, of course, purely superficial. It influences only the rules of drawing, 
i.e. the first item of the diagrammatic rules formulated below. 

Here are the rules for drawing the diagrams and writing down the corresponding analyt- 
ical expressions: 

1. Draw a line corresponding to a propagator of a dot fermion. In Fig. || we used a dashed 
line for its notation. The diagram must include only one such line, since a presence 
of two or more lines of this type implies Qk > 2. It would produce contribution 
to the numerator of equation (|50|) of second or higher order in powers of exp(— /3X). 
The contribution to the current would then be of first or higher order in powers of 
exp(— /3X), which will vanish at A — > oo. Attach m — 1 points to the dot - fermion line, 
where m is the power of Jyk which the contribution of the diagram being considered 
is supposed to have in the expression for the current. These points represent vertices. 
Connect them by lead - electron lines leaving two loose ends corresponding to the 
external operators a\, and au in Eq. floTf). In Fig. |G] we used solid lines to denote 
propagators of lead electrons. In the present work we have considered diagrams for 
m = 2 (Fig. |C) and m = 3 (Fig. |A,B,D). 

2. Introduce the following 2x2 (Keldysh) matrices for every lead - electron and dot - 
fermion line respectively: 



( n n a (a\ ( o CM ^ 



9L(t 



(52) 



The indices "r", "a" and "c" denote the retarded, advanced and correlation Green 
functions which are explicitly defined as follows: 

= (53) 
9 c k Ae) = [l-mtk))\gL(e)-9 a kM 
30 



e (a) M 1 



uj — A ± 27 
= [1 - 2exp(-/3A)] [CM - CM] 

where /(e^) = 1/ [exp(/3efc) + 1] is Fermi function for lead electrons, the factor 
exp(— (3X) is the limiting case of the dot - fermion Fermi function = 1/ [exp(/3A) + 1] 
at large A. As in section |I|the matrices C|52|) emerge as a result of the passage from the 
closed time - path to the normal time path before the transformation to the energy 
representation is carried out. Note that the Lagrange multiplier A appears explicitly 
in the definition of although consideration of energy conservation shows that, in 
fact, it could be absorbed into a shift of u in all Green functions for the dot - fermions. 

3. For every vertex, introduce the following tensor: 



outk in ) 



■I ■/ 



j ko ^js [-^^r 1 ) (-i)X/ » (54) 

where ki n and k out are momenta of incoming and outgoing lead electrons respec- 
tively, and Jk out k in means a matrix element (given by equation (|45|) ). Further, 
Js (W(k out k in )/ (hQ)) are Bessel functions in which the order s is a number of photons 
absorbed in this vertex (clearly, it can be negative), and 



W(k t k- ) 



W, if g L, k in g R 

-W, if k out eR, k in eL ■ (55) 

0, if k ou ti k{ n G L or k ou t, k{ n G R 



It is useful at this point to recall the origin of the appearance of the Bessel functions, 
starting from the expansion of the time - dependent matrix elements Jk utk in {t) (see 
equations fl4~4| ) and (f43|)): 



Jk fk - (t) = J k e^(^ i ») w ^ k i») /Wsill(fit+a) " ,3 ( k »«' fe m )] (56) 
= J k out k m }^ Js\ — ^ — I ( 1) e h \- ( ° vt m > "(*«.t*m)j 
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The symbols (k out ki n ) here, as throughout the present manuscript, denote the leads to 
which ki n and k out belong, that is, LR, RL, LL or RR. The constant phases /3(k out k in ) 
cancel when equation ( |56| ) is substituted into expressions for the diagrams. The tensor 
7] is now a tensor of the fourth rank (unlike in section |I|) because it connects four 
Keldysh matrices. It is written explicitly as, 



Vi 



^ ^ B^i m iBji n i(7 m i n i8 m ' m 8 m i n B m j,B 



"J ■ 



(57) 



m' ,n' ,m,n 



1 1 


' R ~ l = 7S 


d 


-0 


<T <3) - 
> u m'n' 


d 0^ 


. The third Pauli 


v- 1 l ) 


K 1 


1 y 




^0 -1, 





where i? = 



matrix cr^, n , appears here because of the transition from the closed - time path to 
the single - time representation (transformation from integration on the backward 
oriented time branch to integration over normal time axis requires a change of sign). 
Transformation to the physical representation employs the matrices R and Br 1 . 

4. Close around the product of matrices and tensors which has been obtained by the 
following rows (or columns): 2~ 1 / 2 ( — 1, 1) for the g - matrices and 2~ 1 / 2 (1, 1) for the 
£ - matrices (see examples below). The origin of these rows is again the transforma- 
tion from the single - time representation to the physical representation of the Green 
functions and vertices. Multiply the result by the factor, 



Jk'kJn+si + ...+ 



lt...tSm-l 



(58) 



This factor emerges out of the matrix element Jk'kif) in equation (|50|) by means of the 
expansion |56| (while s is replaced by — s). Summation over n is carried out in equation 
(ID below. 



5. Conserve energy in every vertex taking into account absorption of s photons and the 
energy <p dc gained (lost) by an electron going from left to right (from right to left). 

6. Conservation of spin is somewhat delicate. There are two terms in the Hamiltonian 
( f46l) and, correspondingly, two terms in the expression for the current fl50|) . They 
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express spin - flip and normal scattering processes. Therefore, conservation of spin in 
every vertex should take into account two possibilities: 

(a) if incoming spins are opposite, then, spin flip occurs, 

(b) if incoming spins are equal, then, both outgoing spins are the same. 

It is important to stress that the same rules have to be applied to the whole diagram, 
otherwise it can not contribute to the current (|50|) . Namely, 

(a) if incoming spins [p' e and cr'f in Fig. ||) are opposite, then, the outgoing spins 
have to be reversed (i.e. a e = cr'f = —cr' e and 07 = cr' e = —cr'f), 

(b) if incoming spins are equal, then, outgoing spins have to be the same (i.e. cr e = 
cr f = cr' e = a' f ). 

These rules prohibit the spin - flip diagram of type B in Fig. |6|. Indeed, spin - flip in 
the first vertex can happen only if a' e = —cr'f. Then, a ej i = —07,1 = —cr' e = cr'f, and spin 
- flip should occur in the second vertex as well. Consequently, cr e = —Of = cr' e = —cr'f. 
This contribution is therefore forbidden. The only diagram of type B which survives 
is the one with cr e = 07 = o e ^ = 07,1 = a' e = cr' f . On the other hand, four diagrams of 
type A are allowed. Namely, the following equalities a e> i = 07,1, a' e = 07 and cr'f = cr e 
should hold. All the combinations which are left are allowed, that is, a e = ±07 and 

7. Integrate over energy and momentum and sum over spin. Remember that, according 
to Eq. (|50D the current flow starts from the left lead, i.e. k' G L, and ends in the right 
lead, i.e. k e R. The other momenta might belong to any lead: . . . , k m _ 2 G L, R. 
Sum over numbers of absorbed photons si, . . . , s m _i. 

8. Multiply the result by the factor 

_i_ ( _, r -l (i) 2r, ( _ l) F (27r) -( m+ l) (59) 
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where F is the number of closed electron - fermion loops. The first factor in the 
above expression comes from equation (|50|), while the second one reminds us of the 
m — 1 order of the expansion of the S - matrix. Employing Wick's theorem in the 
perturbation series requires pairings of creation - annihilation operators which are then 
expressed in terms of Green functions. This is the origin of the next two factors. The 
last one emerges following Fourier transforms. 

9. Take the limit A — ► oo. 

This completes the list of rules for drawing and calculating diagrams in the time - dependent 
Kondo model. 



C. Explicit expression for the tunneling current 

Let us denote the expression obtained through the above rules for a certain diagram (or 
after summing of expressions for diagrams up to a certain order) as Tr-(ti). The number n 
has the meaning of the total number of photons absorbed or emitted due to the tunneling 
process (it was introduced above in the forth item of the diagrammatic rules). The time - 
dependent current is then given by, 

I +00 

I(t) = -I + cos(nttt + not + arg J n ), (60) 

Z 71=1 

I n = T K (n) + T* K (-n). 

As an example, an expression that is obtained using the above diagrammatic rules for 
diagram A in Fig. (BF) reads, 



1 1 



E Jk'kJkkiJkik'j -p: n 4 / /// de 1 du du' duJt 

fc'GL;fce/?;fcieL,_R;(T,<T',(7i > 



' ^ lim 77TT E (-!> + w' - wi + <f>( kl k) + s 2 hQ) 

\Q«) i',f,i,j,i' v ■ 
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c^i)^<&i)e i V) 



(A 



(61) 



V 1 /,-< 



A useful simplification can be worked out in this equation, noticing that one of the vertices 
connects electron Green functions of the same lead. Indeed, if k\ G L, then k\ and k' belong 
to the same lead, while if k\ G R, then k\ and k belong to the same lead. Therefore, 
according to Eq. ([55]), one has W^ out k in ) = in this vertex. Using the identity, J s (0) = 5 S:0 , 
we can get rid of one of the Bessel functions. Using also the relation J s (— x) = (—l) s J s (x) 
we obtain, 



le 



E 



Jk'kJk^Jkxk' 



lim 



111 



k'£L;k£R;ki£L,R;a,a' ,ai 

dei duj duj' duj\ 



a— (Qk) 4 (2tt) 4 

J2 (-1, l)i0fc$(ei + - + 0( felfc ) + s 2 7ifi) 



(62) 



v 1 / 



v 1 y 



where s = si, if k\ G -R, and s = s^, if &i G L. 

At this stage, the use of a program for analytical calculations (e.g. Mathematical is 
indispensable. The following steps are to be done: (1) the tensor - matrix product is to 
be expanded and expressions for Green functions (|53"D are to be substituted, (2) due to the 
presence of the factor (Qk) in the denominator, only terms of the first order in exp(— (3X) 
must be kept in the numerator (the zero - th order cancel while the terms of higher orders 
are to be omitted), (3) routine integration of the expression obtained through the above 
steps over internal variables e±, u, u' and uo\ is to be carried out. The result is the following: 
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T£(n) 



^ Jk'kJkkiJkik' 
k 1 eL]keR;kieL,R;o-,a' ,cn 

- d> dc - sQ 



(63) 



2^7 



sfl - 17) (e fe - e fcl + (f) (kkl) - s 2 fi - 27) (efc' - e*i + </>(fc'fci) - s x fi + «t) 
/(e*)/( e fa) 



+ 



(efc - ejfe' - dc - sO - 27) (e fc / - e fel + - sj.fi + 27) 

/(e fc ')/( e fci) 



+ 



sfi - fy) ( e fc - e fcl + - s 2 fi - 27 



We assume now that J^k depend only on the leads to which k' and k belong, independently 
of the values of k! and k (hence there are just three coupling - strength constants: Jll, Jrr 
and Jlr-i where Jl(r)r(l) = ^i(-R)^3(i)/l e <i,<j|)- Furthermore, we assume that the densities of 
states in the leads (p(efc), k G L,R), are roughly constant on the energy scale of <fi , fi and 
PU (rigorously speaking, only the combinations p(ek') 1 ^ 2 Jk'k{t)p{^k) \ where k,k' G L, R, 
need to be constant). Then the above expression greatly simplifies since integration over 
one of k, k' or k\ can be carried out. Using also a simple relation, 



Jk'kJkkxJkxk' 



■I 



LR 



2 ~ 



J(kiki)i 



(64) 



we get an expression, 



Ti(n) = -8n\J LR 



PlPr (JllPl + JrrPr 
dek dek 



(65) 



f^k)f(e k ' 



£k — e/c' 



— sfi — ry 



The factor 8 in this equation arises from summation over spins. 

For the diagram B in Fig. |] the above diagrammatic rules, using J s (—x) 
and J s (0) = 8 Si o, yield the following expression: 



-iyj s (x) 



Jk'kJkkiJkik' 



lim 



1 1 1 



dei da; da;' dcui 



a— (Qk) 4 (2vr) 4 

£ l )j9ka^l ~ W + ^1 + 0(fcik) + S 2^) 



(66) 
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Repeating the steps which led from Eq. ( |6^) to Eq. d65|) we obtain, 

PlPr (<7llPl + JrrPr) 
dek dek 



Tf(n) = -27r|j LR 



(67) 



/(efc)/(e*' 



sf2 — ry 



We omitted in this expression some terms which cancel out when substituted into expression 
for the current fl60|) . Note that spin - conservation rules discussed above allow only one 
combination of spins for this diagram, namely, all spin projections are identical. It leads to 
a factor 2 in the above expression instead of 8 in Eq. ( |65f) . 

For the diagram C in Fig. we obtain, using the above diagrammatic rules and J s (—x) = 
(—l) s J s (x) an expression of second order in Jk'kP, 



lim 



1 1 1 



<ie' da; du' 



a-oo (Q^)4(2tt) 3 

£ (-1, l)i<?i£(e' - w + a;' + dc + sftn) 



(68) 



Making the same steps as above we obtain a very simple equation: 
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-{\ 
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1 J 


i' 


i 1 


/ 



T^(n) = -2vr|j LR 



PlP/j 



(69) 



The last term (i.e. i ■ const) cancels when substituted into expression for the current (|BTj). 

In order to calculate contribution of diagram D it is better to stay in the closed time 
- path representation. Then one can easily notice that it can be expressed through the 
contribution of diagram C in the following manner: 
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T £H = T £W / dtp J kiki{t)n kl ; 



(70) 

where n kl = (a^a^J = f(e k ). The matrix elements J kl ki(t) do not depend on the branch 
to which t belongs, i.e. whether i p = t_ or t p — t-, while nfe a does not depend on time at 
all. So the integrals over the normal and the backward oriented time branches cancel. The 
contribution of this diagram to the current simply vanishes. 
Further progress employs the identities 

J s {x)J s+n {x) = 5 nfi (71) 



sJ s {x)J s+n (x) = -(5„,1 + <*n,-l)- 



Substituting Eq. into Eq. ( |S0D and using these equations we obtain for the contribution 
of diagram C to the current (I^ 2 '(t)) the following simple expression: 

j( 2 ) (t) = C 2 [<j) dc + W cos{Qt + a)] , (72) 

2 

where Ci = 1 7r Jl_r Pl/?r. This is the only contribution which is of second order in Jk' k {t)p- 
There appear only DC and the first harmonic. The higher harmonics are not generated in 
this order. 

Finally we sum contributions of all the diagrams and get the following equations for the 
tunneling current: 



I(t) = I (2 \t) + I {3 \t) 
j( 2 ) (t) = C 2 [<p dc + W cos{Qt + a) 



I (3 \~t) = -Io + ^2\I n \cos(nQt + na + a.rgl n 

^ n =l 



In = C3 2-j J* 



fW 



+ Js-n ^)F* 



F(<f>,T,D) 



+00 r+00 

/ de k dt k > 
-d J-d e k - e k 



/(efc)/(efc 



where C- 



2 = 



J, 



LR 



p L p R and C 3 = |107r 



J, 



LR 



(73a) 
(73b) 
(73c) 

(73d) 

(73e) 
(73f) 

PlPr (JllPl + JrrPr), while p L(H ) are 



J s+n (^-)F(<p dc + sQ,T,D) + 



W 



idc 



+ sfl, T, D) 
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densities of states in the leads. The quantities 1^ and 1^ express contributions of second 
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(diagram C) and third (diagrams A and B) orders in Jk'kP respectively. The cutoff D 
is equal to the energy difference between the chemical potential and the bottom of the 
conduction band, while 7 is an infinitesimally small number. Equations (|73|) constitute the 
central formal result of the present section. Equation ( |73t| ) is rather convenient for further 
elaboration (which is our next step) but it can be misleading if a proper care is not taken. 
Indeed, the cutoffs for the two energy integrations are not independent, and their relation 
is to be defined carefully. Consider another expression for F where one of the integrals does 
not require a cutoff: 

m T , d) = -i R e ir ** * m^Mmm±MA _ „t coth » (74 ) 

' ; 2 JJ-D u-e + i-f 2 2 ' v ; 

Here = 1/ (exp[(e — <p)/kT] + 1) and /r(c) = 1/ (exp[e/A;T] + 1) have the meaning of 

Fermi functions in the leads (the left lead being shifted by 0). Some pure imaginary terms, 
linear in 0, that do not contribute to the current have been omitted in this expression. The 
above integral can be written in the same form as in the DC result of Sivan and Wingreen0. 
The physical content of equation (|73d|) for I n is rather transparent: An alternating field 
applied to a two - barrier system actually results in a splitting of the leads energy levels^. 
Therefore the time - dependent current is a result of interference between U DC - like" con- 
tributions, each one of them having an effective bias <p dc + sQ weighted by the appropriate 
product of Bessel functions^. An approximate evaluation of the double integral in the above 
equations is possible both for the linear (0 T) and for the nonlinear (</> ^> T) regimes. It 
yields, 

Re[F(4>, T, *)] = (* ^ * + 26 + ° WkT) + ° (ir/B) l ■' f " << ^ (75, 
I 4> [in ^ + 1 + O {{kT/(j)f) + O {<j)/D)] , if > kT. 

As is well known (see, for instance, Ref. |J), the linear - response perturbation expansion for 
the Kondo problem breaks down at temperatures below a certain value, defined as the Kondo 
temperature Tr-. In our calculation it shows as a logarithmic divergence of the function F 
and, consequently, the amplitude of the current ( |73d[ ) at T — > in the linear response. It 
is remarkable, however, that in the nonlinear case on which we focus our attention here 
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the function F does not diverge with decreasing temperature. Hence, our results for the 
nonlinear response are valid even below the Kondo temperature. This is due to the fact that 
here the non - linear bias plays the role of temperature as the largest low - energy scale. 

When the estimate ([75| ) is employed in equations ([73]) it yields extremely simple expres- 
sions for the current. First, it turns out to be very useful for the static non-equilibrium case 
(W = 0), for which we find, 



C 2 + 0.26C 3 + C 3 In if dc < T, 
^ 2 + C 3 + C 3 ln^, if dc »T. 



(76) 



To the best of our knowledge, such a simple expression for a non - equilibrium tunneling 
current through a Kondo system has not been derived before. We notice that it contains 
the familiar pattern of the zero - bias anomaly (i.e. a peak in the differential conductance 
at zero bias). 

Next, considering expression (|75|) we notice that the large factor \nD appears only in 
a term which is linear in <p. Using equations ( [73 d| ) and ([H]), we then find that this factor 
emerges only in the expressions for the DC and the first harmonic of the AC. It means that, 
in fact, only the direct and the first harmonic are enhanced by the Kondo effect. As for the 
higher harmonics, the interference of the contributions to the current in equations ([TJJ) with 
different effective biases is destructive. This conclusion is further confirmed by numerical 
calculation of the whole spectrum of the tunneling current performed using equations ([73|) . 
In Fig. [TJ the spectra (amplitude of the harmonics I n via their number n, the value of 
1/2 Iq is displayed for the DC) are shown for both the Kondo and the non - interacting 
systems. Values of the parameters used for the calculation are listed in the caption. For the 



non - interacting system we used the equations of Ref. £3] reduced to the case of only one 
resonant level present in the system. Comparison of the top and bottom parts of the figure 
clearly shows a significant enhancement of DC and the first harmonic relative to the others 
in the Kondo system. The ratio of the DC and the first harmonic to the higher harmonics 
might increase even further if D becomes larger. However, the choice of parameters for 
Fig. |7| corresponds to a real quantum dot situation0i. We emphasize that the non - linear 
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response is different from the linear one. Namely, (i) the second and the higher harmonics 
exist although they are small, (ii) the amplitudes of the DC and the first harmonic are not 
determined any more by the ratio D/T as in the linear response, but rather by the ratio of 
D to the largest relevant low - energy scale as we show in the following (see equations ([77D , 
(0) and (§2} below). 

Substituting equation (|75|) into Eq. (|73|), it is possible, in a few limiting cases, to find 
very simple expressions for the DC and the first harmonic of the time - dependent tunneling 
current. Experimentally, the nonlinear Kondo effect is usually revealed as the zero - bias 
anomaly. In a strong alternating field where many - photon absorption is possible (i.e. 
W ^> Q,T), we get for the differential conductance (dl / d(j) ) an expression, 

oi dc „ „ r , ,„ ,„ rt , „ n 



u -c 2 + c 3 ihd/w) + 1] + c 3 — £ p (<p d yn - q , r/n) , if |0 dc | « w, (77) 

oc P zvv \ g \<w/n 

gjdc 

g^- c ^C 2 + Calnp/I^l), if |/ c | » W, 

where P {^(f) dc /Vt — q, T/Oj is the shape function for the side peaks at 4> dc = qQ, q is integer. 
We find that, 

P (V c /fi - q, T/n) « i In (0dc _^ )2 + T2 , if 1/7^ - ?|< 1 and \T/Q\ « 1, while (78) 
P (<p dc /VL - q,T/n) = 0, if \<f) dc /Q - q\ > 1 or |T/fl| > 1. 

In a weak alternating field (W <C f2) , we obtain, 

f)T dc C D 2 /M/\ 2 ' 9 ' 1 

^-ft + f ln WTr5 + C '£y ^f^'-AT,!)), (79) 
where the side - peak shape function P is given by, 

p(^-^,r,c)«^ ( ^_g )a + p - (80) 

Comparing equation ( |77|) with the static expression (76), we notice that the main peak of 
the differential conductance (the one at <p dc = 0) is suppressed in a strong alternating field 
by the factor approximately equal to [ln(D/W) + 1] / [In (D/T) + 0.26]. The side peaks in 
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a strong field are not simple replicas of the central one, as is the case in a weak field. As 
far as the temperature dependence of the zero bias anomaly for the time - independent 
response is concerned, we note that the temperature dependence of the side peaks in a weak 
alternating field is governed by the factor ]n(D/T). In a strong alternating field, however, 
it is determined by the factor In(f2/T). Moreover, the amplitudes of the side peaks in a 
strong alternating field does not decay exponentially with q. In fact, it is roughly constant 
(as long as qVt < W). The half - width of the side peaks in a strong field is approximately 
equal to (fiT) 1 / 2 instead of (DT) 1 ^ 2 . The latter feature is quite favorable for experimental 
observation: since the peaks are rather narrow it is not necessary to go to high frequencies 
in order to resolve them. On the other hand, their magnitude decreases together with the 
ratio Q/T. In Fig. § and Fig. ^| the differential DC- conductance (calculated numerically 
using equations (ff3|)) is shown versus the bias at different magnitudes of the alternating 
field. Suppression of the central peak in the zero bias anomaly is readily manifested in 
both figures. Side peaks do not appear in Fig. ^ where the ratio of the frequency to the 
temperature is not large enough (f2/T = 5). They are well pronounced, however, at a tenfold 
lower temperature (Fig. g). 

Although measurement of an A C with frequencies and amplitudes in the relevant range 
is not an easy task, it might reveal new interesting features of the Kondo effect. Kondo 
contribution to the direct tunneling current is usually revealed in an experiment through a 
special dependence on the parameters (such as In T increase of the conductance or the zero 
- bias anomaly). It was shown above that, as far as the spectrum of the tunneling current 
is concerned, only the DC and the first harmonic are enhanced by the Kondo effect. This 
implies that such kind of parameter dependence can be found in the first harmonic as well 
as in the DC but not in higher harmonics. Employing the approximate expression ( [75]) in 
Eq. ( |73|) we obtain for the first harmonic of the time - dependent tunneling current (I ac ) in 
a strong alternating field (W fi), 

I ac « C 2 W + \c z W [ln(D/W) + 2] , for W > \4> dc \,T, (81a) 



42 



I ac « C 2 W + -C 3 W \ln(D/\(/) dc \) + 1/3] , for W,T < |0 dc |, (81b) 
J ac « C 2 W + C 3 W MD/T) + V 4 ] , f o r W, \<j) dc \ < T. (81c) 

These expressions appear to be very similar to equations ( |77|) for the £>C when |0 dc | and 
are exchanged. Moreover, differentiating J ac with respect to W, we find a peak in the 
differential AC - conductance at zero AC - bias. The shape of the peak is logarithmic at 
large W, while its height is determined by max(\f> dc \,T). We infer that this feature of the 
AC is analogous to the familiar zero bias anomaly in the DC. In the same way that the latter 
is suppressed by the alternating - bias, this "zero alternating - bias anomaly" is suppressed 



by the direct - bias. In Fig. [H] the differential alternating - conductance calculated using 
equations flT3|) is plotted versus alternating - bias W at different values of the direct - bias f> dc . 
Notice the clear similarity with the zero bias anomaly shown in Fig. [8] and its suppression 
with alternating - bias. We do not find any side - peaks in the differential AC- conductance. 
Their traces can be exposed, however, in the dependence of the AC on the direct - bias at 
low temperatures. In Fig. |TT| its derivative (^p?) is drawn versus <p dc for two temperatures. 
At low temperature there appear well - pronounced dips (which are actually peaks in its 
absolute value) at integer multiples of the frequency. 
In a weak alternating field (W <C ft) we find, 

I ac w C 2 W + C 3 W\n{D/\<j) dc \), for T, Q < |0 dc |, (82a) 
I ac w C 2 W + C 3 W [ln(D/0) + 1] , for T, \<j) dc \ < Q , (82b) 
I ac w C 2 W + C 3 W [\n(D/T) + 1/4] , if |0 dc |, ft < T. (82c) 

Inspecting equations (|76|), (|77|) , (^TJ) and (jS^ ) we notice that the values of both direct and 
alternating differential - conductances are basically determined by the logarithm of the ratio 
of D to the largest relevant low - energy scale. 
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V. CONCLUSIONS 



In the present work the problem of non-equilibrium time - dependent electron tunneling 
through an interacting system was studied at some depth. The main attention was focused 
on calculation of time - dependent current in the Kondo regime beyond linear response. A 
tunneling system in this context is naturally described by the time - dependent Anderson 
model. Perturbation expansion of the current within this model, specially adapted for sys- 
tems out of equilibrium was elaborated upon in section O. This formalism combines the 
non-equilibrium Green functions method with a specific approach suggested by Coleman0 
to account for averaging in restricted subspaces which is often encountered in problems 
involving strongly correlated electron systems. We have accomplished the formal part of 
perturbation expansion within this model for the time - dependent case. The task of per- 
forming detailed calculations turns out to be too formidable. Yet, with a slightly elevated 
capability of present day workstations and analytic software programs it should be feasible. 



We then suggested a way to overcome this problem in section |T| where we developed a 
time - dependent version of the Schrieffer - Wolff transformation mapping the time - depen- 
dent Anderson model onto a Kondo - type model. The latter is much easier for treatment 
within perturbation theory. We maintain that it cannot be introduced phenomenologically 
since proper correlation of the time - dependence between the leads and the Kondo coupling 
constant Jkk'if) has to be taken into account. Non-Equilibrium perturbation technique for 
calculation of the tunneling current within the time - dependent Kondo model was worked 



out in section [IV|. Actual calculations were performed up to the third order in J^y (which 
corresponds to sixth order in tunneling matrix elements for the Anderson model) yielding 
extremely simple analytical expressions for the whole spectrum of the tunneling current (see 
equations (|73|), (|76|), fl77|), (|81|) and (|82D). The nonlinear time - dependent current was found 
to be an interference sum of U DC- like" contributions, each one with an effective bias altered 
by the number of absorbed or emitted photons. We stress that our results are valid for the 
non - linear response both below and above the Kondo temperature although for the linear 
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response they are valid only above it. 

There are three novel results in the present research. First, it was found that the Kondo 
effect strongly affects the first harmonic of the alternating tunneling current, no less that it 
affects the DC, while the other harmonics remain relatively small. This result is shown to 
originate from the interference of DC - like contributions to the current (equation ( [75] ) that 
turn out to be rather destructive for all the harmonics except the DC and the first one. The 
higher harmonics are of course generated but their amplitudes are relatively small. This 
result was demonstrated to be remarkably different from that for a non - interacting one - 
level system where all the harmonics emerge together. 

Second, it was found that the zero - bias anomaly in the DC is suppressed by an alter- 
nating field and displays side-peaks at multiples of the basic frequency. This result can be 
easily tested experimentally since it is concerned with measurement of the DC in the Kondo 
regime, which has now been well - established in quantum - dots experiments. Expressions 
([T?|) together with Figs. |8| and || provide an estimate for the preferred range of parameters 
of the system. 

Third, we found a "zero AC - bias anomaly" in the alternating current, i.e. a peak of 
the differential AC - conductance at zero AC - bias. This phenomenon is an AC - analog of 
the familiar zero - bias anomaly of the DC. As in the latter one, it is suppressed by the AC 
- bias, while the former one is suppressed by the DC - bias. There are no side - peaks of 
the differential AC - conductance, however they have a well - pronounced counterpart in the 
derivative of the AC with respect to the DC - bias (in the form of dips). We think that this 
phenomenon shows a pattern of the Kondo effect in the AC yielding a challenging object 
for an experimental search. Equations (|8lD , (^) and Figs. |H], |TT] provide an estimate for 
the necessary range of parameters. Our results on the spectrum of the current indicate that 
effects like this one cannot appear at higher harmonics of the time - dependent tunneling 
current (due to the Kondo effect). We emphasize that expressions (|77|), fl8T|) and (|32|). 
especially the numerical coefficients, are very approximate. They are intended to display the 
basic dependence of the current on the parameters of the system. Quantitative comparison 
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with experiment may be done using the full set of equations (|73|). 

As far as relation to previous relevant works is concerned, we first notice that our ana- 



lytical results for the DC are consistent with the numerical calculations of Ref. |T7[ However, 
being able to consider stronger AC - fields (larger ratio W/Q), we find also an overall sup- 
pression of the zero - bias anomaly, beside the appearance of side peaks. As for the Fourier 
spectrum of the time - dependent tunneling current, we are unable to validate the assump- 
tion suggested in Ref. [17] that all the harmonics beside the DC one can be neglected. Rather, 
the first harmonic is also enhanced, while the second and higher harmonics are generated 
but they are indeed much smaller than the DC and the first oneS Within a specific model, 
some authors^ obtained current spectrum similar to that of a non - interacting system. We 
attribute the difference between this result and ours to a quite peculiar choice of parameters 
used therein. 

We believe that further research on time - dependent aspects of the Kondo effect, in 
particular in quantum dots, is interesting and very timely. Let us mention a few possible 
directions of future research. First, the methods developed in the present work can be 
adapted for solution of the problem of nonlinear response of a Kondo system to a combination 
of alternating magnetic and electric fields. Evolution of the zero - bias anomaly in a magnetic 
field, contrary to its temperature dependence, seems to be the clearest experimentallyfll 
resolved feature of the non-equilibrium Kondo effect. Calculation of complementary effects 
in the time - dependent response and carrying out pertinent experiments seem also to be 
timely. If the effect of a magnetic field can be solely expressed by the Zeeman splitting of the 
energy level in the dot, that is, e a= ±\ = eo ± Ae/2, the technique developed in the present 
work can be easily modified to incorporate it. 

It might also be interesting to consider a non - magnetic (also called "orbital") Kondo 
system^S subject to a strong alternating field. It is believed^! to be realized in some 
recent experiment s0@ in the form of a two - level atomic tunneling system. In this model the 
conducting electrons interact with an impurity atom which can tunnel between two states. 
Tunneling of the atom is assisted by the interaction. At sufficiently low temperatures, the 
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parameters of the model renormalize so that it becomes equivalent to the 2 - channel Kondo 
modelH. Besides the time - dependent shift of the leads which is present also in the 1 - 
channel model and studied in our work, in a two - level atomic system a time - dependent 
field can cause a change of the effective energy separation between the levels. One possible 
effect of such a change is a crossover between Fermi - liquid and non - Fermi - liquid behavior 
(see, for instance, Ref. pq) . It is especially appropriate to point out here that application of 
an AC- field is rather controllable. An appropriate calculation, if followed by an experiment, 
could then further test the hypothesis that the physics of a 2 - channel Kondo model has 
been realized in some experiments,!'!, a point which has been questioned by some authors^. 

Consideration of a multilevel Kondo system in an external alternating field looks very 
attractive. First, it is an appropriate object for quantum dots experiments. Moreover, it 
was arguedEHS that the Kondo temperature in such a system can be enhanced by orders of 
magnitude. This fact could allow an experimental investigation of the strongly correlated 
regime of the Kondo system. On the other hand, application of a time - dependent field to 
a multilevel system leads to highly nonlinear tunneling processes, e.g. resonant frequency 
multiplicationii. Examination of non-equilibrium transport through such a system in the 
Kondo regime might reveal new and interesting effects. 

In closing, we believe that the physics of strongly correlated particles in strong time - 
dependent external fields and restricted geometries looks to be an important and exciting 
subject for further research. 
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FIGURES 

FIG. 1. a) Schematic drawing of the Kondo peaks in the non-equilibrium interacting density of 

states of the dot (impurity). Here is the bare energy level of the dot, while Hur) are chemical 

potentials in the left (right) lead, b) Possible formation of numerous Kondo - resonance peaks in 

the density of states caused by a time - dependent field. The peaks associated with tunneling to 

the left and the right leads are schematically shown near the left and the right barriers respectively. 

$7 is the frequency of the external field. 

FIG. 2. Choice of points on the closed time - path for the external operators in equation (|i3|). 
FIG. 3. Diagrams for the perturbation expansion of the current in the Anderson model. Dashed 

lines stand for slave fermions, dash - dotted lines — for slave bosons, solid lines — for lead electrons. 

A quantity Sj is the number of photons emitted when an electron goes from the dot to a lead, while 

qi is the number of photons absorbed by an electron going from a lead to the dot. 

FIG. 4. Tunneling vertices in the slave - bosons representation of the Anderson model, (a) 

Tunneling of an electron from the dot to a lead is represented as the decay of a slave fermion into 

a slave boson and a lead electron, (b) Reverse process representing tunneling from a lead to the 

dot. 

FIG. 5. Choice of points on the closed time - path for the external operators in equation ([49|). 
FIG. 6. Diagrams for the perturbation expansion of the current in the Kondo model. Solid 

lines stand for lead electrons, dashed lines — for dot electrons. 

FIG. 7. a) Spectrum of the tunneling current for a quantum dot in the Kondo regime (amplitude 

of the harmonics I n via their number n, the value of l/2Io is shown for the dc). The contribution 

of the second order in Jp is not shown. The current is measured in units of C3 ■ kT. W/fl = 4, 

Q/kT = 5, 4> dc /kT = 10, D /kT = 200. b) Spectrum of the tunneling current for a non - interacting 

one - level system with the energy level between the chemical potentials of the left and the right 

leads. The current is measured in units of eT/h, where T is the level width. W/tt = 4, Q/T = 5, 

cf) dc /T = 10, D/T = 200. 

FIG. 8. Differential DC- conductance (in units of C3) versus DC- bias at various values of AC 

- bias W (both are measured in units of kT). 0,/kT = 5, D /kT = 200. 

FIG. 9. The same as in the previous plot but at ten - fold lower temperature. 
FIG. 10. Differential AC - conductance (in units of C3) versus AC - bias at various values of 

DC - bias 4> dc (both are measured in units of kT). Q/kT = 5, D /kT = 200. 
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FIG. 11. Derivative of the AC with respect to DC - bias versus the latter (which is given in 
units of Ml). W/Q. = 2, D/Q = 40. Dashed line is plotted at fi/T = 5, while solid line is at the 
ten - fold lower temperature. 
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